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Abstract 

Taking into account a torsion field gives rise to a negative pressure contribution 
in cosmological dynamics and then to an accelerated behaviour of Hubble fluid. 
The presence of torsion has the same effect of a A-term. We obtain a general exact 
solution which well fits data coming from high redshift supernovae and Sunyaev- 
Zeldovich/X-ray method for the determination of cosmological parameters. On 
the other hand, it is possible to obtain observational constraints on the amount of 
torsion density. A dust dominated Friedmann behaviour is recovered as soon as 
torsion effects are not relevant. 
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1 Introduction 



A generalization of Einstein General Relativity can be obtained by considering a torsion 
tensor different from zero in U4 space-time manifold where connection is not symmetric 
[1, 2]. Such an approach is very interesting today in relation to several extended theories 
of gravity as Superstrings, Supergravity and Kaluza-Klein theories. In particular, torsion 
allows to include spin matter fields in General Relativity and the Einstein-Cartan-Sciama- 
Kibble (ECKS) theory is one of most serious attempt in this direction. However, not all 
the forms of torsion are directly connected to a spin counterpart as it is widely discussed 
in [3]. 

Now if some forms of torsion allow to take into account spin in General Relativity, it 
seems reasonable that they could have had some role into dynamics of the early universe 
when the number density of particles per volume was huge. The presence of torsion 
gives naturally a repulsive contribution to the energy-momentum tensor [4]. In fact, it is 
possible to show that for densities of the order of 10 47 g/cm 3 for electrons and 10 54 g/cm 3 
for protons and neutrons, torsion could give observable consequences if all the spins of the 
particles result aligned. These huge densities can be reached only in the early universe 
so that cosmology is the only viable approach to test torsion effects [4]. However no 
relevant tests confirming the presence of torsion have been found until now and it is 
still an open debate if the space-time is a Riemannian V 4 manifold or not. Considering 
cosmology and, in particular primordial phase transitions and inflation [4, 5, 6], it seems 
very likely that, in some regions of early universe, the presence of local magnetic fields 
could have aligned the spins of particles. At very high densities, this effect could influence 
the evolution of primordial perturbations remaining as an imprint in today observed large 
scale structures. 

^From another point of view, the presence of torsion could give observable effects 
without taking into account clustered matter but resulting as a sort of cosmological con- 
stant. Recent observations seem to point out that the universe is accelerating. Type 
la Supernovae (SNe la) [7], data coming from clusters of galaxies [8] and CMBR in- 
vestigations [9] give observational constraints from which we deduce that the universe 
is spatially flat, low density and dominated by some kind of non-clustered dark energy. 
Such an energy, which is supposed to have dynamics, should be the origin of the observed 
cosmic acceleration. 

In terms of density parameter, we have Q m ~ 0.3,f2A — 0.7, Qk — 0.0, where Q m 
includes non-relativistic baryonic and non-baryonic (dark) matter, A is the dark en- 
ergy (cosmological constant, quintessence,..), k is the curvature parameter of Friedmann- 
Robert son- Walker (FRW) metric. 

Standard matter fluid as source of Einstein-Friedmann cosmological equations gives 
rise to expanding decelerated dynamics. To fit observations, non-standard forms of 
matter-energy have to be taken into account: the net effect should be to implement 
a sort of a cosmological constant which naturally gives rise to a negative pressure capa- 
ble of implementing an accelerated cosmic expansion. On the other hand, Quintessence 
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[10, 11, 12] generalizes this approach taking into account all the mechanisms which give 
rise to negative pressure regimes for cosmic fluid. In particular, scalar fields. 

The cosmological constant problem is one of the main issue of modern physics since its 
value should provide the gravity vacuum state [13], should be connected to the mechanism 
which led the early universe to the today observed large scale structures [14, 15], and 
should predict what will be the fate of the whole universe (no-hair conjecture) [16]. 

In any case, we need a time variation of cosmological constant to get successful in- 
flationary models, to be in agreement with observations, and to obtain a de Sitter stage 
toward the future. In other words, this means that cosmological constant has to acquire 
a great value in early epoch, it has to undergo a phase transition with a graceful exit 
and has to result in a small remnant toward the future coinciding with the observational 
constraints (coincidence problem). 

The today observed accelerated cosmological behaviour should be the result of this 
dynamical process where the present value of cosmological constant is not fixed exactly 
at zero. 

In this context, a fundamental issue is to select the classes of gravitational theories and 
conditions which "naturally" allow to recover an effective cosmological constant without 
considering special initial data. Theories with torsion could match up this point, as we 
will see below, also if the coincidence problem (that is the fine tuning between huge initial 
values of A and very thin today observed constraints) remain essentially unsolved. 

This paper is organized as follows. Sec. 2 is devoted to an essential summary of gravity 
with torsion. In Sec. 3, we show how introducing torsion in cosmological dynamics gives 
rise to a negative pressure extra-term acting as a cosmological constant and we found a 
general cosmological solution. In Sec. 4, we match such a cosmological solution with data 
coming from SNe la surveys and Sunyaev-Zeldovich/X-ray method. Sec. 5 is devoted to 
discussion and conclusions. 



In this section, we give general definitions of torsion and associated quantities. We shall 
use the notation in [1] and in [3]. It is a convention to call U4 a 4-dimensional space- 
time manifold endowed with metric and torsion. The manifolds with metric and without 
torsion are labeled as V 4 . 

Torsion tensor S ab c can be defined, in a U4 space-time manifold, using the antisym- 
metric part of the affine connection coefficients T c ah) that is 



2 Gravity with torsion 





where a, b, c = 0, ... 3. 

In Einstein General Relativity, it is postulated that S ab c = 0. 
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Often in the calculations, torsion occurs in linear combinations as in the contortion 
tensor, defined as 

K a b = —S ah c — S c ab + S b c a — —K a \ , (2) 
and in the modified torsion tensor 

T ab c = S ab c + 25 [a c S b] (3) 

where S a = S ab b . 

According to these definitions, it follows that the afline connection can be written as 

K b = { C ab }-K ab C , (4) 

where { c ab \ are the usual Christoffel symbols of the symmetric connection. 
The torsion contributions to the Riemann tensor can be explicitly given by 

R abc d = R abc d ({}) ~ ^ a K bc d + V b K ac d + KjK b : - K be d KJ (5) 

where R abc d ({}) is the standard Riemann tensor of General Relativity. 

/From Eq.(5), the expressions for the Ricci tensor and the curvature scalar are 

R ab = R ab ({}) - 2V a S c + V b K ac h + KjK bc e - 2S e KJ (6) 

and 

R = R{{}) - AV a S a + K ceb K bce - AS a S a . (7) 

The closest theory to General Relativity containing torsion is the Einstein-Cartan- 
Sciama-Kibble (ECSK) theory. It is described by 

L = V=g(i- + £ m ), (8) 



,2k 

which is the Lagrangian density of General Relativity depending on the metric tensor g ab 
and on the connection T^ b . R is the curvature scalar (7) and C m the Lagrangian density 
of matter fields. 

By the variation of the matter Lagrangian density with respect to the metric we get 

t ab = (9) 
og ab 

which is the symmetric stress-energy tensor while the variation with respect to the con- 
tortion tensor 

is the source of torsion. In many instances, it can be identified with a spin density, but 
there are many cases in which the source of the torsion field is not spin [3]. 
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By the full variation of (8) and introducing the canonical energy-momentum tensor 



S ab = t ab + *y c (r a6c - T bca + T Cab ) , (11) 

where we used the abridged notation *V C := V ' c + 2S cd d , the following field equations are 
derived [1] 

G ab = SnG^ ab , (12) 

and 

T ab c = 8nGr ab c , (13) 

where c = 1, and the symbol V indicates covariant derivative with torsion. 
Eqs.(12) generalize the Einstein equations in a U4. 

Eqs.(13) are algebraic so that it is always possible to cast Eqs.(12) in pure Einstein 
ones, by substituting the torsion terms with their sources. The result is the definition 
of an effective energy-momentum tensor as the source of the Riemannian part of the 
Einstein tensor [1]. In doing so, one obtains 

G ab ({}) = 8nGi a \ (14) 

where G ab ({}) is the Riemannian part of the Einstein tensor. The effective energy- 
momentum tensor is then 



t ab = t ab + 8nG 



a ac bd <-> acd b , cda b 
_4r [d T c] ~ ZT T cd + r T cd 



+ lg ab (K c [d T ed c] + T ecd T ecd ) 

The tensor t ab can be of the form 



(15) 



t ab = (p + p)u a u b -pg ab , (16) 
if standard perfect-fluid matter is considered. 

These considerations are extremely useful for the following, since torsion contributions 
can be treated as additional perfect-fluid terms. 



3 Cosmology with a torsion A-term 

The totally antisymmetric part of torsion can be expressed by a 4-vector 

a a = e abcd S bcd . (17) 

If one imposes to it the symmetries of a background which is homogeneous and isotropic, 
it follows that, in comoving coordinates, only the component o° survives as a function 
depending only on cosmic time (see [3, 17, 18]). 
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For a perfect fluid, the Einstein-Friedmann cosmological equations can be written as 

- = 5~(p + 3p), (18) 

(X o 

and 

/a\ 2 k 8ttG . . 

(a) + ? = — < 19 » 

where a(t) is the cosmological scale factor and k — 0, ±1 is the spatial curvature 
constant. Energy density and pressure can be assumed in the forms [3, 17, 18] 

P = P + f, P = P~f 2 , (20) 

where / is a function related to er°, while p and p are the usual quantities of General 
Relativity. This choice can be pursued since we can define S a j, c = SiaM = f(t) where f(t) 
is a generic function of time which we consider as the source of torsion. For a detailed 
discussion of this point see [3] 

As above, we can define a stress-energy tensor of the form 

tib°* ) = (P + p)u a U b - pQab , (21) 

which, by Eqs. (20), can be splitted as 

j. (tot) . (matter) . . (torsion) /aa\ 

t ab — t ab + t ab ■ 

Due to the contracted Bianchi identity, we have 

tir );b = o ; (23) 
from which, we can assume that (cfr. [19]) 

,(matter);b „ ,(torsion);b „ (OA\ 

t ab — U ■> t ab — U • l^ 4 J 

In the FRW space-time, (23) becomes 

'p + 3H(p + p) = (25) 

which is 

p + 3H(p + p) = -2ff. (26) 
^From Eqs. (24), both sides of (26) vanish independently so that 

fit) — fo — constant. (27) 

In other worlds, a torsion field gives rise to a constant energy density. Taking into account 
standard matter whose equation of state is defined into the above Zeldovich range, we 
obtain 
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Po [-^J + f ■ (28) 

Inserting this result into the cosmological equations, we obtain, in any case, a mono- 
tonic expansion being > 0, d 2 > 0. We obtain an accelerated behaviour if 

p + 3p<2ft, (29) 

so that acceleration depends on the torsion density. 
In a dust-dominated universe, we have 

P = Po(^) 3 + / 2 , P=~fi, (30) 
and then the general solution is 

/ 3 \V3 V3 

a{t) = \i^) [ cosh ^)- 1 ] • ( 31 ) 

Obviously, if fat — > 0, we have cosh(/o£) — (fot) 2 and then a ~ t 2//3 , as it has to be. 

This results tell us that a A-term can be obtained without considering additional 
scalar fields in the dynamics but only assuming that the space-time is U4 (a manifold 
with torsion) instead of V 4 (the standard Riemann manifold of General Relativity). The 
advantage to pass from V4 to U4 is due to the fact that the spin of a particle turns out to 
be related to the torsion just as its mass is responsible for the curvature. From this point 
of view, such a generalization includes the spin fields of matter into the same geometrical 
scheme of General Relativity [1]. 

In what follows, we will compare the obtained solution with data coming from obser- 
vations. 



4 Matching torsion A-term with observations 
4.1 The Supernovae SNe la method 

It is well known that the use of astrophysical standard candles provides a fundamental 
mean of measuring the cosmological parameters. Type la supernovae (SNe la) are the 
best candidates for this aim since they can be accurately calibrated and can be detected 
at enough high red-shift. This fact allows to discriminate among cosmological models. 
To this aim, one can fit a given model to the observed magnitude -redshift relation, 
conveniently expressed as : 

fi(z)=5\og^-d L (z) + 25 (32) 
-"0 



6 



being \i the distance modulus and di(z) the dimensionless luminosity distance. The 
distance in the model we are considering is completely equivalent to the one in a spatially 
flat universe with a non-zero cosmological constant. Thus <1l{z) is simply given as: 

d L (z) = (i + z) f dz'[n M (i + z f + n tm ]- 1 ' 2 . (33) 

Jo 

where Q tor = 1 — Qm plays the same role as the usual Q^. 

The distance modulus, as said, can be obtained from observations of SNela. The 
apparent magnitude m is indeed measured, while the absolute magnitude M may be de- 
duced from template fitting or using the Multi- Color Lightcurve Shape (MLCS) method. 
The distance modulus is then simply /x = M — m. Finally, the redshift z of the super- 
nova can be determined accurately from the host galaxy spectrum or (with a larger 
uncertainty) from the supernova spectrum. 

Our model can be fully characterized by two parameters : the today Hubble constant 
H Q and the matter density Qm- We find their best fit values minimizing the x 2 defined 
as : 

A/ 2/rr o \ \ - ^{ZilHo, Qm) ~ ( o A \ 

X (-Ho, Qm) = 2^ 3 — —-Jl \ 6 V 

i MOi* mz,i 

where the sum is over the data points [21]. In Eq.(34), <r W) is the estimated error of the 
distance modulus and a mz is the dispersion in the distance modulus due to the uncertainty 
a z on the SN redshift. We have: 

5 / 1 dd L \ 

amz = hOO [d-^z-)^ (35) 

where we assume a z = 200 km s^ 1 adding in quadrature 2500 km s _1 for those SNe 
whose redshift is determined from broad features [21]. Note that a mz depends on the 
cosmological parameters so that we will employ an iterative procedure to find the best 
fit values. 

The High-z team and the Supernova Cosmology Project have detected a quite large 
sample of high redshift (z ~ 0.18 — 0.83) SNela, while the Calan-Tololo survey has 
investigated the nearby sources. Using the data in Perlmutter et al. [7] and Riess et al. 
[20], we have compiled a combined sample of 85 SNe as described in detail in [21]. We 
exclude 6 lixely outliers SNe as discussed in [7], thus ending with a total of 79 SNe. The 
results of the fit are presented in Fig.l where we show the 1,2 and 3<r confidence regions 
in the (Qm, Ho) plane. 

The best fit values (with la error) turn out to be : 

Q M = 0.30 ±0.08 , H = 65.4 ± 1.2 km s' 1 Mpc' 1 . 

which allow to conclude that a torsion A-term could explain observation very well. On 
the other hand, we can estimate the torsion density contribution which result to be : 
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fl = (5.6 ± 0.7) x 1(T 30 g cm' 3 , 

which is a good value if compared to the cosmological critical density. It is worthwhile 
to note that, in the case of SNe, the error on Hq does not take into account systematic 
uncertainties due to possible calibration errors. Due to this reason it is so small with 
respect to other results in literature. 

4.2 The Sunyaev-Zeldovich/X-ray method 

Besides the above results, we can discuss how the Hubble constant Hq and the torsion 
density parameter Q tor can be constrained also by the angular diameter distance D A as 
measured using the Sunyaev-Zeldovich effect (SZE) and the thermal bremsstrahlung (X- 
ray brightness data) for galaxy clusters. Distances measurements using SZE and X-ray 
emission from the intracluster medium are based on the fact that these processes de- 
pend on different combinations of some parameter of the clusters (see [22] and references 
therein). The SZE is a result of the inverse Compton scattering of the CMB photons of 
hot electrons of the intracluster gas. The photon number is preserved, but photons gain 
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energy and thus a decrement of the temperature is generated in the Rayleigh- Jeans part 
of the black-body spectrum while an increment rises up in the Wien region. We will limit 
our analysis to the so called thermal or static SZE, which is present in all the clusters, 
neglecting the kinematic effect, which is present in those clusters with a nonzero peculiar 
velocity with respect to the Hubble flow along the line of sight. Typically the thermal 
SZE is an order of magnitude larger than the kinematic one. The shift of temperature 
is: 

AT 



T n 



= y 



x coth ( — 
2 



(36) 



where x = — — is a dimensionless variable, T is the radiation temperature, and y is the 
k B T 

so called Compton parameter, defined as the optical depth r = a T J n e dl times the energy 
gain per scattering: 

y = B ° n e a T dl. (37) 
J m e c 2 

In the Eq. (37), T e is the temperature of the electrons in the intracluster gas, m e is the 
electron mass, n e is the numerical density of the electrons, and <jt is the cross section of 
Thompson electron scattering. We have used the condition T e ^s> T (T e is the order of 
10 7 K and T, which is the CBR temperature is ~ 2.7K). Considering the low frequency 
regime of the Rayleigh- Jeans approximation we obtain 



AT RJ 

Tn 



- -2y (38) 



The next step to quantify the SZE decrement is the one to specify the models for the 
intracluster electron density and temperature distribution. The most commonly used 
model is the so called isothermal (3 model [23]. We have 




M 

2 



(39) 
(40) 

being n eo and T eo , respectively the central electron number density and temperature of 
the intracluster electron gas, r e and (3 are fitting parameters connected with the model 
[24]. For the effect of the cluster modelling see [25]. From (6) we have 

AT _ 2K B a T T eo n eo 

~n~ (41) 

being 

/*oo / / r \ 2\ 2 

dr. (42) 




The integral in Eq. (42) is overestimated since clusters have a finite radius. The effects 
of the finite extension of the cluters are analyzed in [25], [26]. 
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A simple geometrical argument converts the integral in Eq. (42) in angular form, by 
introducing the angular diameter distance, so that 



2 v 1/2-3/3/2 , M _ i 

= = Mi + [fM ^4^— (43) 







2 , 



r(f) 



In terms of the dimensionless angular diameter distances, cIa (such that -D^ = — <1a) 

H 



we get 

AT(fl) _ 2 a T K B T ec n eo 
and, consequently, for the central temperature decrement, we get 




(44) 



AT(fl = 0) 2 q T ^T ec n eo ^ (f |) c 

= -io^^^7(fy^^- (45) 

The factor — d,A in Eq. (45) carrys the dependence on the thermal SZE on both the 
Ho 

cosmological models (through H and the Dyer-Roeder distance (Ia) and the redshift 
(through (Ia)- From Eq. (45), we also note that the central electron number density is 
proportional to the inverse of the angular diameter distance, when calculated through 
SZE measurements. This circumstance allows to determine the distance of cluster, and 
then the Hubble constant, by the measurements of its thermal SZE and its X-ray emission. 

This possibility is based on the different power laws, according to which the decrement 
of the temperature in the SZE, AT (^-°) ; anc i the X-ray emissivity, Sx, scale with respect 
to the electron density. In fact, as above pointed out, the electron density, when calculated 
from SZE data, scales as d~ A x ( nf ZE oc d^ 1 ), while the same one scales as d~ A 2 {nf ~ ray oc 
d^) when calculated from X-ray data. Actually, for the X-ray surface brightness, Sx, 
assuming for the temperature distribution of T e = T e0 , we get the following formula: 

8 * = f^Trhf e ^ dAhx ' (46) 

being 

Jo \n e oJ 

the X-ray structure integral, and ex the spectral emissivity of the gas (which, for T e > 
3xl0 7 , can be approximated by a typical value: ex = , w hh e ~ 3.0xlCT 27 np erg 

cm~ 3 s _1 K^ 1 [24]) . The angular diameter distance can be deduced by eliminating the 
electron density from Eqs. (45) and (46), yielding: 



Sx " e \m e c 2 J eu C H A fl(3/3 - ±, ±, 



(47) 
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where B(a,b) 



T{a)T{b) : 



T(a + b) 
It turns out that 



is the Beta function. 



D, 



^ d ^ (ATp) 2 1 

Hq SxoTeQ @c 



(4* 



where all these quantities are evaluated along the line of sight towards the center of the 
cluster (subscript 0), and 9 C is referred to a characteristic scale of the cluster along the 
line of sight. It is evident that the specific meaning of this scale depends on the density 
model adopted for clusters. In our calculations we are using the so called (3 model. 

Eqs. (47) allows to compute the Hubble constant H , once the redshift of the cluster 
is known and the other cosmological parameters are, in same way, constrained. Since the 
dimensionless Dyer-Roeder distance, oU, depends on Qtor, n m , comparing the estimated 
values with the theoretical formulas for Da, it is possible to obtain information about 
Q m Qtor, and H . Recently distances of 18 clusters with redshift ranging from z ~ 0.14 
to z ~ 0.78 have been determined from a likelihood joint analysis of SZE and X-ray 
observation (see Table 7 [27] and reference therein). Modeling the intracluster g clS clS db 
spherical isothermal /3-model allows to obtain constraints on the Hubble constant Hq in a 
standard A-FRW model. We perform a similar analysis using angular diameter distances 
measurements for a sample of 44 clusters, constituted by the 18 above quoted clusters 
and other 24 already know data (see [22]). In Fig (2), our data are plotted against two 
theoretical models. 



h=0 . 7 



Q 




Figure 2: We plot the measure angular diameter distances with the error bars and two 
theoretical model: a Einstein- De Sitter model (bottom curve) and a flat A-FRW model 
with fl m = 0.3, Q\ = 0.7 (upper curve), both with Hq = 70 km s _1 Mpc -1 



As indicated in [22, 27], the errors a are only of statistical nature. Taking into account 
our model with torsion (31), the theoretical expression for the angular diameter distances 
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D A is 

D A (z) = j-^—d^z) = [ Z dz'[Q M (l + z'f + fi^]- 1 / 2 

(i + zy 1 + ,2 jo 

We find the best fit values for Q tor and ifo, minimizing the reduced x 2 '- 



X 2 (^0,^m)=E 



a 



(49) 



(50) 



The results are shown in the Fig. (3), where we see the contours corresponding to the 
68.5% and 98% confidence levels: The best fit values (at la) turn out to be : 



90 



80 



X 70 



60 



50 




0.2 0.4 0.6 0.8 1 

Figure 3: 1, 68.5% and 98% confidence regions in the (flu, H ) plane. The best fit values 
are : Vt M = 0.3, H = 68 km s" 1 Mpc -1 . 



tt M = 0.30 ±0.3 , H = 68 ± 6 km s~ l Mpc~ l 
in good agreement with the above fit derived from SNe la data. 
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5 Discussion and conclusions 

In this paper, we dealt with a cosmological model where torsion is present into dynamics. 
After the discussion of how such a contribution modifies cosmological Friedmann-Einstein 
equations, we have shown that the net effect of torsion is the introduction of an extra-term 
into fluid matter density and pressure which is capable of giving rise to an accelerated 
behaviour of the cosmic fluid. Being such a term a constant, we can consider it a sort 
of torsion A-term. If the standard fluid matter is dust, we can exactly solve dynamics 
which is in agreement with the usual Friedmann model (to be precise Einstein-de Sitter) 
as soon as torsion contribution approaches to zero. 

The next step is to compare the result with observations in order to see if such a 
torsion cosmology gives rise to a coherent picture. We have used SNe la data, Sunyaev- 
Zeldovich effect and X-ray emission from galaxy clusters. Using our model, we are capable 
to reproduce the best fit values of H and Qm which gives a cosmological model dominated 
by a cosmological A-term. In other words, it seems that introducing torsion (and then 
spins) in dynamics allows to explain in a natural way the presence of cosmological constant 
or a generic form of dark energy without the introduction of exotic scalar fields. Besides, 
as we have seen in Sec. 4, observations allows to estimate torsion density which can be 
comparable to other forms of matter energy (~ 5.5 x 10~ 30 g cm~ 3 ). 

However, we have to say that we used only a particular form of torsion and the 
argument can be more general if extended to all the forms of torsion [3]. Furthermore, 
being in our case the torsion contribution a constant density, it is not possible to solve 
coincidence and fine tuning problems. To address these issues we need a form of torsion 
evolving with time. This will be the topic of a forthcoming paper. 
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